Activation and protonation of dinitrogen at the FeMo cofactor of nitrogenase 
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The protonation of N2 bound to the active center of nitrogenase has been investigated using state-of-the-art 
density-functional theory calculations. Dinitrogen in the bridging mode is activated by forming two bonds to Fe 
sites, which results in a reduction of the energy for the first hydrogen transfer by 123 kJ/mol. The axial binding 
mode with open sulfur bridge is less reactive by 30 kJ/mol and the energetic ordering of the axial and bridged 
binding modes is reversed in favor of the bridging dinitrogen during the first protonation. Protonation of the 
central ligand is thermodynamically favorable but kinetically hindered. If the central ligand is protonated, the 
proton is transferred to dinitrogen following the second protonation. Protonation of dinitrogen at the Mo site 
does not lead to low-energy intermediates. 



I. INTRODUCTION 

Although dinitrogen is the main part of the atmosphere, 
it is, in its molecular form, inaccessible to biological organ- 
isms. Its inactivity is caused by the triple bond — one of the 
strongest covalent bonds in nature. While high pressure and 
high temperature are required to convert N2 to NH3 in the 
industrial Haber-Bosch process, biological nitrogen fixation 
achieves the same goal at ambient conditions. For this pur- 
pose nature employs the enzyme nitrogenase, which is one of 
the most complex bioinorganic catalysts in nature. 

Nitrogenase converts N2 to biologically accessible 
ammonia. 1,2 * 3 * 4 During the reaction, non- stoichiometric 5 
amounts of hydrogen are produced. 

N 2 + (6 + 2x)H+ + (6 + 2x)e~ -> 2NH 3 + xU 2 . 

Values for x range from <1 (Ref. 6) upwards and there is 
an ongoing discussion about the question whether hydrogen 
is produced in stoichiometric amounts. 7 In the gas phase, the 
first protonation is by far the most difficult step in the reac- 
tion. The catalyst has to find a way to dramatically reduce this 
barrier. 

The active center of the enzyme, shown in Fig. 1, is the 
FeMo cofactor, MoFerSgN-homocitrate. The FeMo cofactor 
is linked to the protein via two amino acid residues. The ques- 
tion arises why nature employs such a large multi-center clus- 
ter. 

The reaction mechanism has being studied for over 40 
years, but the atomistic mechanism of substrate conversion at 
the FeMo cofactor still remains an open issue. In a previous 
paper, we presented the most salient features of the reaction 
cycle as it emerged from our simulations. 8 In the present pa- 
per, we present the results of state-of-the-art first principles 
calculations of the most difficult step in the reaction pathway, 
the first protonation of N2 bound to the FeMo cofactor. 

Nitrogenase consists of two proteins: (1) the molybdenum- 
iron protein, which holds the FeMo cofactor (FeMoco) and 
(2) the iron protein, which hydrolyzes MgATP and uses the 
energy obtained to provide the molybdenum-iron protein with 
electrons. 
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FIG. 1 : The FeMo cofactor of nitrogenase with its homocitrate ligand 
and the two residues linking the cofactor to the protein: histidine and 
cysteine. 



The structures of the two proteins were resolved already 
in 1992. 9 * 10,11 Since then increasingly more refined crystallo- 
graphic data became availabl e) 12 - 13,14 until, in 2002, a central 
ligand of the FeMoco, shown in Fig. 1 was found, 15 which was 
undetected in previous studies. While crystallographic studies 
determined the central ligand to be either a C, N, or O atom, 
theoretical studies iyil&AS provide strong support for it to be 
a nitrogen atom. It turns out that the central ligand (N x ) plays 
a critical role in the mechanism^ 

Kinetic studies of the mechanism of biological nitrogen 
fi xa ti on 20 j 2i i 22 i 23 i n( ii ca t e that the rate-limiting step of the re- 
action is the dissociation of the two proteins. In each of 
these association-dissociation cycles one electron is trans- 
ferred to FeMoco. The so-called Thorneley-Lowe scheme 
provides insight into the first reduction steps, stating that N2 
binds after three or four electrons have been transferred to the 
MoFe protein. Theoretical models 24,25 indicate that geometri- 
cal changes of the backbone of the Fe-protein are responsible 
for using the energy from MgATP to transfer electrons to the 
MoFe protein. 

Since the appearance of the first crystal structures, 
nitrogenase has been subject to numerous theoretical 
investigations, 2 ^ 2 ^ 

47 i 48 i 49 i 50 i 5i i 52 Unfortunately, the central ligand was un- 
known at that time, so that the conclusions have to be 
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reconsidered. However, several more recent theoretical 
studieAiSii2ii2ii2iSi5i have been performed taking the central 
ligand into account. 

A critical parameter for theoretical calculations is the ox- 
idation state of the cofactor. By comparing theoretical 
results 1718 * 54 with experimental observations, a consensus has 
been reached that the cofactor in the resting state is charge 
neutral, i.e., [MoFe 7 S 9 N] . 

The mechanisms for nitrogen fixation proposed up to now 
may be divided in two main groups: (1) conversion at Mo and 
(2) conversion at Fe. 

• The Mo atom has been the target of numerous 
theoreticatfdSiiSi&& studies. A strong indication in fa- 
vor of Mo comes from experiment: an Mo-based model 
complex 57 has been found to catalytically reduce N2. 
However, there are also active nitrogenases 3 with Mo 
replaced by V or Fe, which indicates that Mo is not es- 
sential. 

• Numerous proposals have been made for the reduction 
involving the Fe atoms. (1) The direct way for N2 to 
bind at Fe atoms is head-on binding. 30 - 31 - 38 *S 4 AS In this 
mechanism N2 binds to one Fe atom of the central clus- 
ter, where it is protonated until one and then the sec- 
ond ammonia dissociate. The intermediates of the com- 
plete cycle have been identified, albeit still without the 
central ligand by Rod and co-workers^iSiffi Recently, 
Hinnemann and N0rsko , \iS extended this proposal to 
the cofactor with central ligand. The theoretical pre- 
dictions are, however, at variance with this mechanism, 
because the calculations show that N2 binding is highly 
endothermic. 58 (2) Sellmann and co-workers 59 * 60 - 61 - 62 
suggested an opening of the cage in analogy to smaller 
Fe complexes. In his model, two octahedrally coordi- 
nated low-spin Fe atoms positioned in close proximity 
bind dinitrogen between them, where it is reduced. The 
mechanism has been lined out up to the first four pro- 
tonations. (3) Our own recent calculations^ support 
the view of cage opening even though, in our model, 
the Fe atoms are in a high- spin tetrahedral coordina- 
tion, which points towards a quite different chemistry 
than expected for low- spin octahedral Fe atoms. The 
opening and closing of an SH bridge at the cofactor is 
critical for the activation of N2 and for the dissociation 
of the second ammonia. A critical role is attributed to 
the central ligand, which enables required bond rear- 
rangements. In a similar approach, we were able to 
explain many experimental findings in the conversion 
of acetylene by nitrogenase. 63 (4) Another proposal 19 
suggests the opening of a sulfur bridge upon coordina- 
tion of water to an Fe atom, complete protonation of 
the central ligand, and dissociation of ammonia. Then 
N2 inserts into the central cavity of the cofactor, where 
one nitrogen atom is fully protonated and dissociated, 
which closes the catalytic cycle. This intriguing pro- 
posal seems to be in conflict with isotope exchange (ES- 
EEM/ENDOR) experiments 64 that exclude an exchange 
of a central nitrogen ligand. 



II. CALCULATIONAL DETAILS 

We considered the complete FeMo cofactor with truncated 
ligands as in the previous study. 54 The histidine was replaced 
by imidazole, the homocitrate ligand by glycolate and the cys- 
teine, bound to the terminal iron atom by an SH group. 

We performed density-functional theoryS^ (DFT) calcu- 
lations based on the projector- augmented wav o 67 i 68 (PAW) 
method. The gradient-corrected PBE (Ref. 69) functional 
was used for exchange and correlation. The PAW method is a 
frozen-core all-electron method. Like other plane- wave-based 
methods, the PAW method leads to the occurrence of artifi- 
cial periodic images of the structures. This effect was avoided 
by explicit subtraction of the electrostatic interaction between 
them. 70 Wave-function overlap was avoided by choosing the 
unit cell large enough to keep a distance of more than 6 A 
between atoms belonging to different periodic images. We 
used a plane wave cutoff of 30 Ry for the auxiliary wave func- 
tions of the PAW method. The following shells were treated 
in the frozen-core approximation: Fe [Ne], Mo [Ar3d 10 ], S 
[Ne], O [He], N [He], and C [He]. The following sets of pro- 
jector functions were employed, Fe 2s2p2d, Mo 2s2p2d, S 
2s2p2d, O 2s2pld, N 2s2pld, C 2s2pld, and H 2slp, which 
provides the number of projector functions per angular mo- 
mentum magnetic and spin quantum numbers m, s in each 
main angular momentum channel L 

Atomic structures were optimized by damped Car- 
Parrinellc 71 molecular dynamics with all degrees of freedom 
relaxed. The convergence was tested by monitoring if the ki- 
netic temperature remains below 5 K during a simulation of 
0.05 ps (200 time steps). During that simulation, no friction 
was applied to the atomic motion and the friction on the wave 
function dynamics was chosen sufficiently low to avoid a no- 
ticeable effect on the atomic motion. 

Transition states were determined by applying a one- 
dimensional constraint on the atomic positions. In the present 
application, bond-length, angle, and torsion constraints were 
used. To get a first upper bound for the barrier, the specific 
constraint was varied within 1000 molecular-dynamics (MD) 
steps. If this upper bound is less than 20kJ/mol, the barrier 
will be easily overcome and it has not been calculated more 
accurately. In case of a higher estimate, the bond length was 
fixed to discrete values around the transition state to maximize 
the energy while all unconstrained degrees of freedom were 
allowed to relax to minimize the energy. Proof that this ap- 
proach, when converged, determines exactly first-order tran- 
sition states is given elsewhere^ 

The reaction rates Y can be estimated using Y = 
To e~ EA ^ kuT ^ from the calculated activation energy Ea and 
a typical vibrational frequency To = 3 x 10 13 s _1 correspond- 
ing to about 1000 cm -1 . 

The FeMoco exhibits seven high- spin iron atoms anti- 
ferromagnetically coupled to each other. Many different 
spin configurations may easily lead to metastable states in 
conventional collinear spin-polarized calculations. There- 
fore, we used a noncollinear description of the spin den- 
sity for our calculations. In a noncollinear description each 
one-electron wave function is a two-component spinor wave 



3 



function This method not only correctly describes 

real noncollinear spin states, which occur within the reac- 
tion mechanism, but also avoids artificial barriers between 
different spin configurations occurring in collinear calcula- 
tions. Our resulting spin distribution is therefore independent 
of the random starting conditions. Such a dependence is a 
common problem of conventional (collinear) spin-polarized 
calculations for this system. We found that the spin order- 
ing depends on subtle changes in the atomic structure. It also 
changes between different states of the reaction mechanism. 
The spin orderings encountered in our calculations are given 
in Fig. 2, where we follow the notation introduced by Lovell 
et al4£ noncollinear spin arrangements have been found, in 
this study, only for energetically unfavorable states, which is 
why we do not specify them further. 

The spin quantum number S is specified alongside with the 
corresponding structures in Figs. 3 and 4. A spin with S = 1, 
for example, corresponds to a triplet. Noncollinear states are 
indicated by a spin quantum number that differs from half- 
integer values. This corresponds, in analogy with the collinear 
expression, to the absolute value of the integrated spin density 
divided by H. 




BS6 BS7 



FIG. 2: The two relevant spin orderings BS6 and BS7 as well as the 
spin ordering of BO. BS6 and BS7 differ by a spin flip of Fe3. BS7 is 
the same spin ordering as obtained for the cofactor in the resting state 
M. It is characterized by ferromagnetically coupled Fe-Fe and Fe- 
Mo pairs, which are antiferromagnetically aligned relative to each 
other (Ref. 54). The spin ordering only defines the directions of the 
site spins, not their absolute value. 

During the reaction, protons and electrons are transferred 
to the cofactor and the substrate. In this work we made the as- 
sumption that electrons and protons are transferred in a corre- 
lated way, i.e., that one proton is transferred with each electron 
transfer. This assumption implies one of the two scenarios: 
either a reduction of the cofactor increases the proton affinity 
such that a proton transfer is induced, or, if the proton trans- 
fer precedes the electron transfer, it implies that the electron 
affinity is sufficiently enhanced by the positive charge next to 
the cofactor to induce an electron transfer to the cofactor. This 
is the main assumption in our work besides the accuracy of the 
density functionals and the neglect of the protein environment. 
This assumption has been shown to be valid for the cofactor 
before the binding of the substrate^A 

The energies of protons and electrons, which are consumed 
during the reaction, affect the overall reaction energy. Thus 
we need to define a value that reflects their energies in the en- 



vironment. For protons, the relevant environment is the proton 
transfer channel, while for electrons it is expected to be the P- 
cluster. The exact energies cannot be determined by theory 
alone. 

As a consequence of our assumption that the reduction and 
protonation occur in a correlated manner, only the sum /in of 
the energies of protons and electrons is relevant for the rela- 
tive energies of the intermediates. A range of possible values 
can be derived by comparing experimental x-ray and extended 
x-ray-absorption fine structure (EXAFS) data with our calcu- 
lated geometries: we found indirect evidence for the cofactor 
being unprotonated in the resting state, while being protonated 
in the reduced stated Therefore, /in is sufficiently high to 
drive protonation, that is, /i H > ^(MH) - E(M). On the 
other hand, no protonation occurs under the same conditions, 
but in the absence of MgATP. Thus the chemical potential in 
the absence of MgATP, denoted by /i^, must be sufficiently 
low not to drive protonation, that is, /ijj < E(MR) — E(M). 
As two MgATP are hydrolyzed in each electron transfer, the 
difference between the chemical potentials with and with- 
out MgATP is smaller than twice the energy of hydrolysis of 
MgATP, that is, /i H - < 64.4kJ/moli2i It is smaller be- 
cause a fraction of the energy supplied by MgATP will be dis- 
sipated. Therefore, we assume the lower bound for /in, that 
is, /in = ^(MH) — E(M), in our calculations. This is the 
most conservative assumption possible. A less conservative 
value would make the reactions including protonation more 
exothermic. 

Previous studies 19,38,39,40,58 made a different choice for /in, 
namely, /in = ^(fb). This would be the appropriate 
choice if the protons and electrons would be obtained from 
gaseous hydrogen. While the production of gaseous hydrogen 
2H + + 2e~ -^H 2 is energetically neutral with this choice, this 
reaction is exothermic by 7 1 kJ/mol when using our choice of 
/iH = ^E(H.2 ) + 35 kJ/mol. Thus our reaction energies can be 
translated to the values for H2 reference by adding 35 kJ/mol 
per added hydrogen atom. We additionally list the reaction 
energies with H2 as the reference energy in parentheses after 
the values we obtain with our /in for all those energies were 
influenced by the choice of /in- 

In our work we evaluate not only the energetics of the inter- 
mediates, but also the barriers for the transitions. To estimate 
the barriers for protonation, we need, however, to specify a 
proton donor, which models the proton channel. We used the 
ammonium ion to mimic the proton donor. Note, however, 
that this choice only affects the barriers but not the relative 
energies of the intermediates. Our finding that the barriers for 
protonations are small and will therefore be overcome easily 
indicates that the barriers for protonation are not relevant for 
the overall picture. 

The notation for the structures is chosen as follows. The 
complexes of dinitrogen with the cofactors are given in letters 
in alphabetic order according to the number of proton transfers 
and in numerals for their energetic order. A numeral denotes 
the ground state for the selected composition. 
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III. RESULTS 



A. Protonation of the central ligand 




A1 

2.5 BS6 



FIG. 3: N2 binding modes, their spin states, and their spin ordering 
(BS state) obtained from DFT. 

The most difficult step of the reaction from N2 to NH 3 is the 
first protonation of dinitrogen. In the reaction in the gas phase 
we find that the reaction step N2 + H + + e~ — > N2H is en- 
dothermic with 164kJ/mol using the choice for /in described 
under calculational details (and 199kJ/mol with H2 as refer- 
ence). The main goal for nitrogenase must be to lower this 
barrier. At the FeMo cofactor, the reaction energy is dramati- 
cally reduced to only 41 kJ/mol (76kJ/mol). Nevertheless, in 
the catalyzed reaction this remains the most endothermic step. 

We explored eight different isomers of the N2H adduct 
bound to the MH 2 state of the cofactor. MH 2 represents the 
resting state (M) of the cofactor reduced by two electrons and 
protonated at two of the sulfur bridges. The isomers are shown 
in Fig. 4 and will be discussed in the following. 




BO 

S = 4 
-37 kJ/mol 



B1 

S = 2 BS7 
kJ/mol 



B2 
S = 1 BS6 
19 kJ/mol 



B3 
S = 2 BS7 
24 kJ/mol 




B4 

S = 2 BS7 
29 kJ/mol 



B5 

S = 1.029 
32 kJ/mol 



B6 

S = 0.873 
48 kJ/mol 



B7 
S = 1 BS6 
108 kJ/mol 



FIG. 4: Different binding modes of the N2H adduct at the FeMoco. 
The spin state and the spin ordering (BS state) are given, as well as 
the energies, which are given relative to Bl. 



Nitrogen in the bridged binding mode A I opens the struc- 
ture of the cofactor and leaves space for the coordination at the 
central ligand. Protonation of the central ligand leads to the 
most stable isomer with one protonated nitrogen, namely BO, 
which is lower in energy by 37 kJ/mol than the second most 
stable configuration, the protonated dinitrogen bridge Bl. 

There is only limited space for a donor to access the cen- 
tral ligand in Al. Correspondingly, we find a large barrier 
for protonation of 79kJ/mol from NH4 , which was used as 
the model for the proton donor. This barrier corresponds to 
a reaction rate of the order of 0.1-1 s _1 , which is substan- 
tially lower than the electron transfer rate. Furthermore, as 
discussed below, Al, the common starting configuration for 
both BO and Bl can be protonated with a negligible barrier at 
dinitrogen leading to Bl. Since the pathways to BO and Bl 
are in direct competition, protonation of the central ligand is 
kinetically hindered. Therefore, we consider a reaction path 
via protonation of the central ligand unlikely. 

In BO, one iron atom is in a planar threefold coordination 
shell, which appears to be an unusual configuration. To ex- 
plore if this structure is an artifact of our structural model, we 
investigated if this iron atom could form a complex bond to 
a nearby water molecule, which would restore a tetrahedral 
coordination of the Fe atom. However, in our calculations 
including an additional water molecule, no significant com- 
plex bonds between the three-coordinate Fe atom and a water 
are formed: for BO an Fe-0 distance of 2.7 A was calculated, 
which is substantially larger than typical complex bonds to Fe. 

One may ask if the solvent effects affect the barrier to ac- 
cess BO. We find that the interaction between an additional 
H2O molecule and the triangular Fe atom in the initial state 
A I and in the transition state is even weaker than in the fi- 
nal state BO, expressed in even longer Fe-0 distances. From 
the absence of strong interactions we conclude that the large 
barrier is not appreciably influenced by the solvent. 




FIG. 5: Intermediates for the protonation of the central ligand. 

A slow rate of formation does not rule out that the cofactor 
is accidentally trapped in this low-energy intermediate. There- 
fore, we investigated the subsequent steps starting from BO. 
We find that after one reduction the next proton attaches to the 
bridged dinitrogen, which is similar to the second most stable 
intermediate discussed below. However, upon protonation of 
dinitrogen, the proton at the central cluster is destabilized and 
is transferred to the unprotonated nitrogen atom of dinitrogen. 
The reaction energy for the internal proton transfer from the 



5 



TABLE I: Bond distances (A), angles (°), and vibration numbers 
(cm -1 ) of the N2 binding modes and free N2. They are compared 
to experimental data on N2 and a model complex. 





AO 


theory 
Al 


N 2 


N 2 


experiment 
Model Complex" 


d(N-N) 


1.146 


1.173 


1.105 


1.098 


1.182(5) 


d(N-Fe7) 


1.777 


1.827 






1.770(5) 


d(N-Fe3) 




1.834 






1.779(5) 


d(SH-Fe3) 


2.298 


2.349 








Z(N-N-Fe7) 


176 


147 






180 


Z(N-N-Fe3) 




147 






180 


P(N-N) 


1979 


1792 


2339 


2331 


1778 



a Reference 1781 



central ligand to the protonated dinitrogen is —69 kJ/mol. The 
reaction has a barrier of 20 kJ/mol. Thus even if the central 
ligand is protonated, the reaction mechanism quickly leads 
back to what we consider the most likely pathway. Interme- 
diate structures of this side branch of the reaction cycle are 
shown in Fig. 5. 

B. Protonation of bridged dinitrogen 

The second most stable isomer with one protonated nitro- 
gen is the bridged binding mode Bl shown in Fig. 4. 

As shown previously, 54 dinitrogen can dock at the FeMo co- 
factor MH2 in two configurations with similar energies. They 
are shown in Fig. 3. The axial binding mode, denoted by AO, 
is slightly more stable than the bridged one, denoted by Al. 

One of the most interesting questions of biological nitrogen 
fixation is how dinitrogen is activated for the first protonation. 
Therefore, we investigated the two binding modes in some 
detail. 

The most relevant geometric parameters and the N-N 
stretching vibration frequencies of dinitrogen bound to the co- 
factor are given in Table [0 and compared to the experimen- 
tal data obtained from a model complex with a N2 bridging 
two low-coordinated iron atoms. 78 The elongation of the N-N 
bond as well as the reduction of the vibration frequency with 
respect to gaseous N2 are indications of the activation of the 
triple bond for the following protonation. The elongation of 
the dinitrogen bond and the reduction of the stretch frequency 
in A I compares well with the model complex, while AO ap- 
pears to be much less activated. 

The activation of A I can be traced back to occupied N2- 
7r* orbitals as shown in Fig. 6. These orbitals can be seen as 
an antisymmetric combination of the Fe-N bonds. Interest- 
ingly they are only dominant in the minority spin direction 
of the two neighboring Fe sites Fe3 and Fe7. The origin is 
that the interaction with the unoccupied d states split the 7r* 
orbital in a bonding and an antibonding orbital. The bond- 
ing orbital, still having 7r* character but containing the Fe-N 
bonding contribution, becomes occupied and is located about 
1 eV below the highest occupied molecular orbital (HOMO), 



while the antibonding orbital, having the Fe-N antibonding 
contribution, remains unoccupied. The former, bonding or- 
bital is the relevant frontier orbital for the protonation. 




1 J 1 1 1 u 1 I lLi 1 1 

-8 -6 -4 -2 2 4 



Enenrgy (eV) 




Enenrgy (eV) 



FIG. 6: Density of states of the bridged (Al, top) and the axial (AO, 
bottom) binding modes illustrating N2 activation. The full line in- 
dicates the density of states projected onto the d orbitals of the Fe 
atoms directly bound to N2 . The shaded curve is the density of states 
projected onto the tv and 7r* orbitals of dinitrogen. The insets show 
the relevant wave functions. Note that both states have occupied 7r* 
orbitals, which is more dominant and lower in energy in the bridged 
binding mode A J. 

The energy to add a hydrogen atom to dinitrogen 
in the bridged binding mode A I and to obtain Bl is 
41kJ/mol (76kJ/mol), substantially less than the 164kJ/mol 
(199kJ/mol) of the gas-phase reaction. is energetically 
slightly more favorable by 8 kJ/mol than protonating the other 
nitrogen atom of N2 . 

The protonation of the reduced complex from an ammo- 
nium, which mimics the proton donor in our calculations, pro- 
ceeds with a negligible barrier of only 4 kJ/mol and is exother- 
mic by 63 kJ/mol. Note, however, that the calculated protona- 
tion energies taken individually are not as reliable as the reac- 
tion energy, as the former depend on the choice of NH4 as the 
proton donor. Nevertheless, this calculation indicates that pro- 
tonation of the reduced complex at dinitrogen proceeds with- 
out difficulty and much more readily than protonation at the 
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central ligand. 

Interestingly, the spin ordering of the cofactor, as obtained 
from DFT, changed from BS6 in Al to BS7 in Bl during 
the reduction and protonation. (See Fig. 2 for a definition of 
the spin labels.) This corresponds to a spin flip of the iron 
atom Fe3, to which N2 as well as the SH group are bound. 
In comparison, the FeMo cofactor in the resting state has the 
same spin arrangement BS7 as BL After two protonations 
and reductions, BS7 of the resting state is transformed to a 
noncollinear spin distribution, which changes to BS7 during 
dinitrogen binding in the axial mode AO and to BS6 dur- 
ing the conversion to the bridged binding mode Al. These 
transitions show the importance of allowing spin flips and 
noncollinear spins when simulating the reaction. Previous 
calculations^&iSdiy 8 assumed a fixed spin ordering during the 
entire reaction, which may cause errors in the energy profile. 

C. Protonation of axial dinitrogen 

If N2 is axially bound as in AO during the first protona- 
tion, the lowest-energy protonation site is the terminal nitro- 
gen atom. The resulting structure is B2, shown in Fig. 4. The 
protonation does not induce major structural changes. The 
energy of B2 is 19kJ/mol higher than the energy of B 1, i.e., 
the bridged mode. Nevertheless, even though the energy for 
protonation of the reduced complex from an ammonium is, 
compared to the one of the bridged mode, 30kJ/mol smaller, 
the reaction is still exothermic with 33kJ/mol and proceeds 
with a negligible barrier of less than 10 kJ/mol. 

Also in this case, the spin ordering changed during the re- 
duction and protonation from BS7 in AO to BS6 in B2. 

D. Other binding modes 

Only 5 kJ/mol above the energy of E2, we find the complex 
B3 bridging two Fe atoms with a single nitrogen atom. This 
structural principle is found again later in the reaction cycle. 8 
However, like the axial binding mode, it lies substantially, that 
is, by 24 kJ/mol, higher in energy than the bridged mode BL 

The protonated dinitrogen bound to molybdenum, i.e., B4, 
lies 29kJ/mol above BL Together with B7, it will be dis- 
cussed in Sec. IIII El where we discuss the potential role of 
molybdenum. 

As in B5, dinitrogen can also bridge the two Fe atoms with 
its axis perpendicular to the Fe-Fe direction, so that both ni- 
trogen atoms are connected to both Fe atoms. A similar bind- 
ing mode has recently been found experimentally for nitrogen 
bridging two zirconium centers. 79 The energy of this structure, 
i.e., B5, lies 32kJ/mol above that of BL 

In the relevant intermediate B6 of the mechanism proposed 
by Hinnemann and N0rskov, 58 dinitrogen binds axially to one 
Fe atom like our structure B2. In contrast to B2, however, the 
sulfur bridge is still intact and the bond between this Fe atom 
and the central ligand is broken. This structure is 48 kJ/mol 
above Bl and it is 29 kJ/mol less stable than the axial binding 
mode B2 discussed previously. 



E. Molybdenum as coordination site 

Numerous theoretical^i 4 ^ 4 ^ 4 ^ 4 ^^ and experimental 57,80 
studies have been performed to investigate the Mo atom as a 
possible adsorption site. For this reason, we discuss nitrogen 
coordination to molybdenum here in some detail. As shown 
previously, 54 the molybdenum atom is, according to our sim- 
ulations, not a favorable adsorption site for nitrogen: nitrogen 
bound to molybdenum is higher in energy by 50 kJ/mol than 
the one bound to the iron atoms. Nevertheless we investigate 
the most favorable pathway for the first protonation via dini- 
trogen bound to the Mo site. 

As shown previously, nitrogen adsorption at pentacoordi- 
nate Mo is endothermic by 30kJ/mol. 54 This indicates that, 
even if the coordination site is vacant, dinitrogen binds to Mo 
only for fleetingly short periods of time. If protonation pro- 
ceeds sufficiently easy so that the proton is transferred during 
these short periods, dinitrogen may be stabilized bound to Mo. 
However, protonation leading to B4, shown in Fig. 4 is ener- 
getically unfavorable. B4 is 29 kJ/mol higher in energy than 
Bl, with N2H bound to Fe atoms. 

Durrant 45 proposed a transition of the protonated dinitro- 
gen B4 into a bridging position between Mo and Fe as in B 7 
shown in Fig. 4. We find B 7 to be 79 kJ/mol higher in energy 
than B4. 



IV. DISCUSSION 

The following picture emerges from our calculations: dini- 
trogen exists in two competing docking modes at the cofac- 
tor, the axial mode AO and the bridged mode A 1. 54 They 
are separated by a large barrier of 66kJ/mol, which, how- 
ever, is still sufficiently small to be overcome. 54 Before the re- 
duction and protonation take place, the axially bound dinitro- 
gen is even slightly more stable by 6 kJ/mol than the bridged 
configuration. 54 We assume that axial and bridged modes are 
in equilibrium until the proton is transferred. However, the 
electron transfer and the protonation reverse the energetic or- 
der between them so that the axially bound dinitrogen B2 ends 
up 19 kJ/mol higher in energy than the bridged dinitrogen BL 

In the bridged binding mode Al, N2 is activated for accept- 
ing a proton by forming bonds to the Fe sites. As a result the 
7r* orbitals are occupied, which implies that the triple bond 
is effectively broken. The occupied 7r* orbital exposes two 
nucleophilic lobes to which a proton can bind. 

These occupied 7r* orbitals are binding combinations of 
the dinitrogen 7r* orbital with the d orbitals of the Fe sites. 
They are only occupied in the minority spin direction of the 
Fe atoms binding to dinitrogen. This is because only the high- 
lying d orbitals of the minority spin direction can lead to a sta- 
bilization of the 7r* orbital; hybridization with the low-lying 
Fe-d states of the majority spin direction would shift the 7r* 
orbitals further up, while the corresponding bonding orbital is 
mostly of Fe character. Note here that both Fe sites binding to 
dinitrogen have the same spin orientation in AL 

In the axial binding mode AO, the corresponding 7r* orbital 
lies close to the Fermi level and its weight on dinitrogen is 
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substantially smaller than in the corresponding orbital of the 
bridging mode. These differences between the bridged and 
axial binding modes are also reflected in the bond-length ex- 
pansion and the reduction of the stretch frequency of the dini- 
trogen bond, and suggest a smaller activation of the axial bind- 
ing mode. 

Nevertheless, dinitrogen is readily protonated both in the 
axial and in the bridged binding mode after the transfer of one 
electron. Compared to the gas phase, the reaction energy for 
the first protonation is dramatically reduced from 164kJ/mol 
(199kJ/mol) for the gas phase to 41 kJ/mol (76kJ/mol) at the 
FeMo cof actor in the bridged mode. 

Interestingly, the most stable site for protonation is the cen- 
tral ligand. The resulting isomer is more stable by 37 kJ/mol 
than the one with protonated dinitrogen. However, this iso- 
mer is not considered relevant because protonation of the cen- 
tral ligand directly competes with protonation of dinitrogen. 
While protonation of the central ligand requires a large bar- 
rier of 79kJ/mol to be overcome, protonation of dinitrogen 
proceeds nearly barrier-less. We attribute the large barrier to 
the breaking of the bond between the central ligand and one 
of the Fe atoms, which changes its tetrahedral coordination to 
a trigonal planar one. 

Even if the central ligand is protonated, the second proton 
is added to dinitrogen and induces an internal proton transfer 
from the central ligand to dinitrogen. 8 Thus even if the first 
proton binds to the central ligand, this step neither poisons the 
catalyst nor does it lead to an entirely different branch of the 
reaction cycle. This side branch may potentially be relevant at 
small turnover frequencies. 

The branch via protonation of the central ligand is reminis- 
cent of the proposal by Huniar et al. p 9 who propose complete 
protonation of the central ligand and its cleavage as ammo- 
nia, before dinitrogen binds to the cofactor. In their study, 
the central cage is opened by a water molecule coordinating 
to one Fe atom. Additional calculations would be necessary 
to directly compare our results with their proposal. However, 
the latter seems to be in conflict with isotope exchange (ES- 
EEM/ENDOR) experiments 64 that exclude an exchange of a 
central nitrogen ligand. 

A mechanism via dinitrogen bound to the Mo site is in- 
consistent with our calculations. Binding to the Mo atom is 
substantially more endothermic than binding to the Fe sites. 

Dinitrogen bridging two Fe atoms is part of the proposal 
by Sellmann and co-workers 59 * 60,61,62 . Our model differs from 
theirs, in that all Fe atoms remain tetrahedrally coordinated 
and in the high- spin state. We are not aware of any theoretical 
investigations of this system including the cofactor with its 
central ligand. 

The presence of the central ligand is crucial for the reac- 
tion mechanism of nitrogen fixation at the FeMo cofactor, and 
may explain why the mechanism has remained elusive for a 
long time. The main feature of the central ligand is its ability 
to form a variable number of bonds to the six Fe atoms. The 
central ligand changes its coordination from sixfold to fivefold 
and fourfold. This allows other ligands such as nitrogen and 
sulfur to form and cleave bonds to the Fe atoms without de- 
viations from the preferred tetrahedral coordination of the lat- 



ter. Tetrahedral coordination seems to be a common structural 
principle to nearly all relevant intermediates of the reaction 
cycle. This is particularly apparent when dinitrogen binds: 54 
while the Fe atom to which dinitrogen binds in the axial bind- 
ing mode AO maintains its coordination shell by giving up its 
sulfur bridge, the bridged binding mode Al would result in an 
unfavorable fivefold coordination of one Fe atom, if the latter 
would not give up its bond to the central ligand. 

Let us mention here some variants of the reaction steps dis- 
cussed: our calculations do not allow to distinguish between 
Fe3 and Fe7 as potential docking sites of dinitrogen. We find 
that the axial binding mode AO is more stable, that is, by 
30kJ/mol, on Fe7 than on Fe3. The bridged mode accessed 
via the axial binding at Fe7 is more stable by 14kJ/mol than 
the other variant. We expect that these differences are sensi- 
tive to the environment of the cofactor: due to the large mo- 
tion of the sulfur atom, the opening of the SH bridge may be 
strongly affected by the shape and the specific interactions of 
the cavity. Thus, no conclusive answer regarding the initial 
binding site can be given at this point. 

The cofactor has an approximate threefold symmetry, 
which is broken by the ligands and the protein environment. 
As long as the protein environment is not taken into account, 
as in the present study, the energetics will proceed similar for 
all three equivalent orientations. Nevertheless, the position of 
the proton channel in the protein indicates that nitrogen fixa- 
tion proceeds near the iron atoms Fe3 and Fe7. Furthermore 
the cavity in this region provides sufficient space to accom- 
modate nitrogen bound to an Fe atom. 54 

V. SUMMARY 

We have studied the first protonation of dinitrogen at the 
cofactor on the basis of density-functional calculations. We 
made an effort to explore the phase space for the reaction with- 
out prejudice for one particular model of the mechanism. A 
large number of intermediates and the barriers between them 
have been explored and placed into perspective. While an un- 
ambiguous determination of this reaction step is not yet possi- 
ble, three possible branches could be identified. One proceeds 
via dinitrogen in the bridged configuration, a second one pro- 
ceeds via axially bound dinitrogen, and the third proceeds via 
protonation of the central ligand. The latter, however, is un- 
likely to play a role in the reaction cycle due to kinetic compe- 
tition with direct protonation of the bridging dinitrogen. The 
activation of dinitrogen is discussed in detail. 
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